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ABSTRACT 



Context. Direct imaging of circumstellar disks requires high-contrast and high-resolution techniques. The angular differential imaging 
(hereafter ADI) technique is one of them, initially developed for point-like sources but now increasingly applied to extended objects 
such as disks. This new field of application raises many questions because the disk images reduced with ADI depend strongly on the 
amplitude of field rotation and the ADI data reduction strategy. Both of them directly affect the disk observable properties. 
Aims. Our aim is to characterize the applicability and biases of some ADI data reduction strategies for different disk morphologies. 
A particular emphasis is placed on parameters mostly used for disks such as their surface brightness distribution, their width if the 
disk is a ring, and local features such as gaps or asymmetries. We first present a general method for predicting and quantifying those 
biases. In a second step we illustrate them for some widely used ADI algorithms applied to typical debris disk morphologies: inclined 
rings with various inner/outer slopes and width. Last, our aim is also to propose improvements of classical ADI to limit the biases on 
extended objects. 

Methods. Simulated fake disks seen under various observing conditions were used to reduce ADI data and quantify the resulting 
biases. These conclusions are complemented by previous results from NaCo L' real-disk images of HR4796A. 
Results. As expected, ADI induces flux losses on disks. This makes this technique appropriate only for low- to medium-inclination 
disks. A theoretical criterion is derived to predict the amount of flux loss for a given disk morphology, and quantitative estimates of 
the biases are given in some specific configurations. These biases alter the disk observable properties, such as the slopes of the disk 
surface brightness or the radial/azimuthal extent of the disk. Additionally, this work demonstrates that ADI can very easily create 
artificial features without involving astrophysical processes. For example, a particularly striking feature appears for a ring when the 
amplitude of field rotation is too small. The two ring ansae are surrounded by two flux-depleted regions, which makes them appear 
as bright blobs. This observation does not require any astrophysical process such as dust blown by radiation pressure, as proposed 
previously in H-band images of HR 4796A. 

Conclusions. The ADI techniques behave as spatial filtering algorithms and can bias disk observables. Therefore, the filtering process 
needs to be properly calibrated when deriving disk parameters from processed images. 

Key words, high angular resolution - data analysis - circumstellar matter 



1. Introduction 

The study of circumstellar disks is essential for understand- 
ing the formation of planetary systems. Direct imaging of 
these disks has revealed asymmetries, warps, gaps, troncatures, 
density waves, or other features that are the results of in- 
teractions between the disk and its environment. About 160 
disks are now resolved from their visible to thermal emission 
(http://www.circumstellardisks.org, Stapelfeldt). Debris disks 
are a particularly interesting class of optically thin disks since 
planets are already formed, if any, and faint structures in the dust 
di stribution can betray the ir presence, as proposed for HR 4796A 
in lLagrange et all d2012bl) . 

The vicinity of the star and its lumin osity however Urnit 
the performance of both ground-based dRacine et al] [l999) 
and space-based dSchneider& Silverstond u003l) high-contrast 
imaging because of bright quasi-static speckles. Slowly evolv- 
ing, they add up over time and eventually become the dominant 
noise source at separa tions below a few arcs econds, depending 
on the star brightness dMacintosh et al.ll2005l) 



The principle of differential imaging is to subtract a refer- 
ence f rame from the targ et image to reduce quasi-static speckle 
noise (M arois et alj|2006h . This reference fra me is sometimes 
also c alled a reference Point Spread Function (lLafreniere et al.l 
2007), hereafter PSF, but we will use here the more generic term 
reference frame. This frame can be obtained in various ways: ei- 
ther with a reference star or with the target itself observed at a 
different field of view orientation (ADI), a different wavelength 
(spectral differential imaging), or polarization (polarimetric dif- 
ferential imaging). This paper focuses on ADI. 



Angular differential ima ging has alrea d y ach ieved 
it results on debris disks: Buenzli et al. (2010) fori] 



signifi- 
cant results on debris disks: IBuenzli et al] (2010) for in stance 
detail ed the morpholog y of the Moth. iBoccaletti et al.l d2012l) 
and ICurrie et a .1 (120 12) reve aled the inner part of HD 32297, 
iThalmann et al. d201 ll) and Lagrange et al. d2012bl) studied the 
ring of HR4796A. lLagrange et al.l d2012c ) revealed the inner 
part of the /3 Pictoris disk and constrained the project ed planet 
position relative to the disk, and lThalmann et alj (2010) resolved 
the gap in the transitional disk around LkCa 15. However, side 
effects might alter the apparent morphology of the disk and 
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its observable properties: flux self-subtraction, change in disks 
slo pes, width, etc. The model-matching procedure described 
by iBoccaletti et ail d2012f) to capture the true morphology of 
HD 32297 is a good illustration of how difficult it is to dismiss 
ADI artifacts. 

Angular differential imaging uses the differential rotation 
between the field of view and the pupil occurring on an alt- 
azimuthal mount to distinguish between the speckle halo and 
any on-sky source. When the pupil tracking mode is used, the 
field rotates at the same rate as the parallactic angle while the 
pupil is stable. The parallactic angle is the angle between the 
great circle through the object and the zenith, and the hour circle 
of the object. 

When extended objects such as disks are imaged, the chal- 
lenge of ADI is to build reference frames with a speckle pat- 
tern highly correlated to the target image but without capturing 
the flux of the disk. If the disk flux is captured in the reference 
frame, the reference subtraction will decrease the disk flux in the 
residual image: this is the problem of self-subtraction. Not only 
does it lead to a lower signal-to-noise ratio (S/N) for the disk but 
it also biases the observable parameters of the disk. 

The objective of this study is twofold: 

1 . To guide the observing strategy and data reduction by pro- 
viding key figures to know how much flux loss is expected 
for a given disk geometry and field of view rotation. A theo- 
retical criterion is derived to quantify this loss for the specific 
case of an edge-on disk. 

2. To highlight the observable parameters that can be trusted 
from ADI-based disk images and point out the artifacts po- 
tentially created by this technique. 

We will first describe the simulation procedure, and in 
Section 3 will build theoretical criteria to analyze the applica- 
bility of ADI to disks. A qualitative (Section|4|i and quantitative 
(Section|5]) description of the biases induced by ADI is then pre- 
sented before reviewing specific ADI algorithms more adapted 
for disk reductions in Section[6j Systematically investigating the 
parameter space of these algorithms to derive the most accurate 
disk parameters is beyond the scope of that section. 

2. ADI algorithms and simulation procedure 

Several data reduction strategies to treat ADI data have been de- 
veloped over the past few years. They differ in the way the im- 
ages are selected and combined to build a reference frame. The 
selection can use all available images or only a reduced sample 
based on a disk binary mask, depending on how much the par- 
allactic angle has changed from the working frame, how close 
in time the image is from the working frame, or how correlated 
they are. The combination can be a simple median, an average, 
or a linear combination of all selected frames. 

We focus here on three specific algorithms referred to as 
classic ADI (cADI), ra dial API (rADI) and L OCI. These pro- 
cedures are described in Lagrange et al.l d2012af) . We will review 
three additional techniques specific for disks that use iterations, 
binary masks for the disk (mcADI, mLOCI), and damped LOCI 
(dLOCI). 

- In cADI (Marois et al the reference frame is a simple 
median or mean of all individual images of the cube. The key 
parameter is the total amount of field rotation available, AO, 
typically 20° to 80° for observations of one to several hours. 

- In rADI, the reference frame is e stimated radially in annuli as 
presented in M arois et alJ d2006l) . The key parameters are the 



separation criterion Ng, usually 1 to 5 FWHM, the number of 
images N used to build one reference frame (between 5 and 
30) and the radial extent of the annulus dr. For each working 
frame, the reference frame is obtained by taking the mean 
of the selected N images. The selection of the N images is 
made temporally, after sorting the images according to their 
hour angle and selecting the N closest to the current working 
frame for which the field has rotated more than the separation 
criterion Ng. 

- In LOCI dLafreniere et al.ll2.Q07h . the reference frame is es- 
timated locally in subtraction areas using linear combina- 
tions of all sufficiently rotated frames. The linear coeffi- 
cients, called c , are computed to minimize the residuals in 
an optimization area. The geometry of the optimization area 
is an annulus section parametrized through the aspect ratio 
g, and the surface expressed in number of PSF cores A^. The 
geometry of the subtraction area is parametrized through the 
annuli width dr. One key parameter is the separation crite- 
rion N& as in rADI. 

- McADI and mLOCI are variants of cADI and LOCI using 
masks to limit the contamination of the reference frame by 
the disk. They are detailed in Section [53] 

- DLOCI is an algorithm detailed in Pue vo et all (|2012). It is 
applied here to disks and is discussed in Section lo31 

To study the artifacts generated by those algorithms, we gen- 
erated fake disk images, inserted them into data cubes at the 
appropriate parallactic angle and convolved them with a syn- 
thetic or real telescope PSF. In some cases, noise coming from 
real data was also added to these simulation data cubes. The 
fake disks were assumed to be axisymmetric and were generated 
from an initial dust volumetric d ensity p(r,z) = ppR( r)Z(r,z) 
with the simulation code GRaTer ( Augerea u et al.lfl999h . R(r) = 

{/ r \-2au i \-2q-„„,)- - 5 

\\t) + \k) [ i s tne ra di a l dust distribution in the 
mid-plane and Z(r, z) is the vertical dust distribution at radius 
r. 



3. Applicability of ADI to disks: the problem of 
self-subtraction 

When the disk azimuthal extent is large with respect to the field 
rotation, the reference frame used to subtract the light of the cen- 
tral star and reveal circumstellar matter contains some flux from 
the disk and self-subtraction occurs. 

3.1. Pole-on disks 

An axisymmetric disk seen pole-on is entirely self-subtracted, 
because the surface brightness distribution (hereafter SBD) is 
symmetrical with respect to the central star. In this configu- 
ration, ADI is unsuitable and the PSF subtraction should be 
made using a reference star observed close to the target, prefer- 
ably a t the same parallac tic angle, to maximize speckle corre- 
lation dBeuzit et al.lll997h . Nevertheless, non-axisymmetric fea- 
tures such as bright blobs can still be revealed through ADI, but 
this application remains limited. 

3.2. Edge-on disks 

The situation is different for disks seen edge-on such as the /? 
Pictoris or AU Mic disk. The key parameter is the vertical extent 
of the disk at a separation x relative to the field rotation at that 
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Final position 
of the disk in 
.the data cube 



A6/2 / 




Initial position 
of the disk in ' 
the data cube 



Fig. 1. Uniform disk of height 2ho seen edge-on at the beginning, 
middle and end of the observing sequence. In red, the final disk 
image after cADI processing, in blue the part of the disk that is 
entirely self-subtracted. 



radius. In the very theoretical case of a uniform disk of vertical 
scale height ho, with a total field rotation AO, and a reference 
frame built by median-combining all available images (cADI), 
the residual disk is entirely self-subtracted for a separation x < 

~ Mn(Afl/4)' as iU ustrate d m Figure Q] 

For a non-uniform disk, with a vertical profile parametrized 
with an exponential index y so that I(z) = e , there is no 
longer a sharp boundary radius Xu m between a totally preserved 
region and a totally self-subtracted region. However, we can 
identify three different regimes, depending of the scale height 
h at the given radius with respect to AO. These are illustrated 
in Figure [2] for an initial disk of vertical scale height h, with an 
exponential index y — 2 (Gaussian vertical profile) and a field 
rotation AO = 60°. Either the disk is well preserved in its whole 
vertical extent (for a radius beyond 7 x h) , or the vertical extent 
is reduced while the mid-plane brightness is well preserved (be- 
tween 7 and 6 x h) or the mid-plane is also dimmer (below 6 x h). 
The last two regimes are illustrated by vertical profiles shown in 
linear and magnitude scales in the last four graphs of Figure [2] 
These graphs show that brightness profiles can vary significantly 
whether the reference is combined by median or mean. With a 
mean, the flux in the mid-plane is higher (top chart of left graph 
of Figure |2]i and the vertical extent of the disk is also larger but 
negative regions appear on each side of the mid-plane, as visible 
in linear scale in the middle left graph of Figure [2] A median 
prevents negative regions from appearing on noiseless data, or 
limits their extension in the presence of noise and is more robust 
to noise when many frames are available for reference frame 
combination. Unless otherwise specified, the reference frame is 
built by median combination in cADI throughout. 

A theoretical criterion to quantify the disk flux loss can 
be derived in this specific vertical brightness distribution. Let 
Aa = arctan be the angular vertical height of the disk seen 
from the star. We aim to characterize the remaining disk flux 
at the point of vertical height z and radius r, corresponding 
to an angular height i = arctanyj^j. The remaining disk flux 
/ at a given position in the reduced cADI image, normalized 
to the initial disk flux (for a reference frame built by median- 
combination), follows the exponential law, 



f( Kc ) = \-e K < with k c (i) 



(A6/4y 



(Aa)? 



(1) 



For instance, in the disk mid-plane (i = 0°), k c = ( <A f ^ ) T and 
for k c = 3, 95% of the disk flux is preserved. The derivation of 
that law is detailed in AppendixlAland is verified experimentally. 





Fig. 2. Top: radial intensity profile along the disk mid-plane of 
the initial disk (black) and cADI images (reference frame built 
by median-combination in red and by mean-combination in pur- 
ple). Middle: two vertical intensity profiles at the radii corre- 
sponding to the green dashed lines of the top figure, in linear 
scale. Bottom: same vertical profiles in magnitude scale. 

A similar criterion can be derived for an edge-on disk with 
exponential vertical profile processed with rADI. In this case, 
however, the field rotation term A6 is replaced by a term de- 
pending on the separation criterion N$ and the number of images 
N used to build the reference PSF If we call Ae = x N 
the mean field rotation corresponding to N images and <x = 
arctan ( n s xfwhm ^ ^ se p arat j on criterion expressed in terms of 
angle, the disk flux / in the mid-plane follows the law 



f(K r ) — 1 - e with K r 



3.3. General case 



a- + Ae/2 
Aa 



(2) 



In the two previous subsections, we showed that ADI is inappro- 
priate for pole-on disks and we derived the key criterion k that 
allows one to predict the amount of flux losses for an axisym- 
metric edge-on disk. 

In the general case of a non-axisymmetric disk with a ran- 
dom inclination i, the amount of flux loss cannot be summa- 
rized by this straight-forward mathematical criterion. Instead, 
we show below that ADI can be seen as a spatial frequency fil- 
tering algorithm, and we use this study to answer the question to 
what extent ADI can be applied. 

ADI, a spatial frequency filtering algorithm: To show that cADI 
can be seen as a filtering algorithm, we built a frequency filter 
based on Fourier analysis and applied it to disk models to com- 
pare the results with cADI-processed images. The tool uses as 
input a simulated disk image and an amplitude of field rotation 
AO and performs the following steps for each pixel of coordinate 
(r, 0) in the image: 

- it extracts the lD-azimuthal profile at the radius r 

- it sets to zero the points of the profile whose distance to is 
greater than AO/2 

- it builds the Fourier spectrum of the resulting profile and re- 
moves from this spectrum any frequencies lower than a pre- 
defined cutoff frequency uq 
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- last, it reconstructs the azimuthal intensity profile by an in- 
verse Fourier transform of the spectrum, and reads the inten- 
sity at the pixel position (r, 6). 

Different values for the cutoff frequency uq were tested and the 
choice of ^ leads to an image close to the cADI processed im- 
age, when the reference frame is estimated by mean combina- 
tion. This means that cADI can be seen as a low-pass filter of 
the azimuthal frequencies of an image, the cutoff frequency be- 
ing that corresponding to the total amplitude of field rotation A6. 

A comparison between disk images processed by cADI and 
by our frequency-filtering algorithm is given in Figure [3] for a 
theoretical model of the debris disk HD 141569. We used three 
different amplitudes of field rotations: 30°, 50° and 80°. 

The disk inclination is 52°, which is the critical value below 
which ADI induces so much self-subtraction to make it inappro- 
priate. The scattered light model of that disk was made using the 
GRaTer code, starting from a radial distribution for the disk opti- 
cal depth, using a disk ellipticity of 0. 1, an anisotropic scattering 
factor of -0.2, a Gaussian vertical profile, and a constant open- 
ing angle. The initial model is displayed at the top of Figure [3] 
The disk contains several rings and anisotropic scattering makes 
the side inclined toward the line of sight brighter. The image of 
the disk processed by cADI shows that the morphology of the 
two rings has been altered. They are dimmer and the brighter re- 
gion of those two rings is thinner, distorted and not any longer 
the brightest region. Increasing the total amount of field rotation 
yields more disk flux but does not prevent the distortions cre- 
ated in the disks. All these observations are also valid for the 
images created by the filtering algorithm. Therefore we can con- 
clude that they result from filtering the low spatial frequencies 
that have an orthoradial component. Areas with negative flux be- 
tween the two rings also appear. These negative areas are better 
visible in the brightness azimuthal profile at r — 162AU shown 
in dashed lines in the plot of Figure|4] They are well reproduced 
by the filtering algorithm (plain lines in Figure @J, confirming 
the validity of the interpretation of cADI in terms of frequency 
filter. 



4. Qualitative and morphological effects of ADI on 
disks 

4.1. Negative flux artifacts: extinctions, blobs and 
asymmetries 

Most of the artifacts created in ADI are due to over-subtracted 
areas that can make fake structures appear. 

To highlight those effects, we reduced simulated images of 
the disk of HR4796 as if they were observed in the same con- 
ditions as in the real ob servations on April, 7 2010, detailed in 
Lagr ange et al.l (1201 2bl) . The same field of view rotation (77.5°) 
and the same number of images in the datacube (103) after sort- 
ing poor-quality images were used. The evolution of the paral- 
lactic angle for the 103 images of the data cube is plotted in 
Figure [5] at the top. 

We used five additional parallactic angle evolution patterns: 

- the same parallactic angles as previously but with a disconti- 
nuity corresponding to eliminating 21 images from the data 
cubes at three different positions within the cube (see Figure 
[5]at the top.) This discontinuity can mimic data sorting due to 
poor atmospheric conditions or quality of the adaptive optic 
correction 
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W Frequency filtering 
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Fig. 3. HD 141569 disk model (top), and a comparison between 
cADI processed images with different field rotations and resid- 
ual flux as predicted by the frequency filter. The color scale is 
identical and linear for all images. 



- a reduced field of view rotation of 23° observed after cul- 
mination of HR4796 in the sky, as in Thalm ann et al.l (1201 lh 
(see Figure|5]at the bottom.) 

- a reduced field of view rotation of 23° observed before cul- 
mination 

The initial disk is convolved with a synthetic PSF in L' . The 
reductions were performed with cADI, rADI and LOCI (with a 
separation criterion N s of 1.5 FWHM and 3 FWHM). 
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1.0 




Angle in radians 

Fig. 4. Comparison between azimuthal brightness profiles of 
HD141569 at a radius of 162 AU, initially and after various 
cADI processing or frequency filtering. rad corresponds to the 
right panel in Figure [3] 



4.1 .1 . General effects of cADI 

The first striking effect is the extinction along the minor axis of 
the disk seen in Figure [6] At this distance, the amplitude of field 
rotation is low and self-subtraction is significant, resulting in a 
progressive decrease in the width of the annulus as we approach 
the minor axis. 

With a 77.5° field rotation, no other artificial feature appears, 
except for a slightly negative region inside the annulus and two 
point-like sources along the minor axis with a flux 91% lower 
than the ansae. Nothing distinguishes this image from those with 
discontinuous parallactic angle evolutions: the introduced breaks 
have undetectable effects in cADlQ. 

For a reduced field rotation, two prominent artifacts appear 
for the northeast and southwest ansae of the annulus (left and 
right side in Figure [6] respectively). There are two extinctions 
on each side of the ansae, resulting in two blobs on the parts 
of the annulus that are most distant from the center. Applying 
the frequency-tool described in the previous section shows in- 
deed that the azimuthal spatial frequencies at the location of 
the extinctions are of the same order of magnitude as the cut- 
off frequency due to ADI-filtering, which causes the enhanced 
flux losses in that region. 

As mentioned in Lagra nge et al.l d2012bl) . there is no need 
to involve radiation pressure-induced dust blow-out to explain 
the appearance of thes e bright structures named streamers by 
iThalmann et al.l d201 lb in their high S/N data along the semi- 
major axis. The radial profiles along the semi-major axis at the 
bottom of Figure|6]is unambiguous: the initial outer edge profile 
(small black dashed line) is well reproduced after cADI reduc- 
tion (green line) even with a reduced parallactic angle amplitude 
(brown line) along the major-axis Q. 



1 down to a detection level of 0. 1 % with respect to the ansae intensity 

2 A reduced parallactic angle amplitude attenuates the outer edge in- 
tensity by a uniform factor of 18% but does not alter the shape and log- 
arithmic slope of the profile along the semi-major axis. Only the inner 
edge is altered because of the presence of a negative area with a min- 
imum intensity at 0.95". This conclusion does not hold any more in a 
direction 11° away from the semi-major axis: in that case even the outer 
edge radial profile is altered, as shown by the violet curve in Figure|6] 
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Fig. 5. Top: parallactic angle evolution of HR4796 as observed 
on April, 7 2010. The dashed lines show the artificial discon- 
tinuities introduced for test purposes. Bottom: parallactic angle 
evolution for a total field of view rotation A9 = 23° before and 
after the star culmination in the sky. 

A deprojected view of the disk is helpful to understand where 
the c ADI azimuthal profile deviates from the original disk profile 
(Figure |7). The extinctions described quantitatively below are 
shown here by the local minimum in the azimuthal profiles (plain 
and dotted red lines). They are not present with a 77.5° field 
rotation (green lines). 

4.1 .2. Effects specific to rADI 

In rADI (Figure |SJ, the annulus appears thinner than in cADI. 
The blobs are still present at the tip of the ansae. The negative 
regions are emphasized. The differences between the two sets of 
parallactic angle evolution, observed and reduced amplitude, are 
not as high as in cADI, because the criterion that rules the refer- 
ence frame selection is now Ng rather than AG. We see two neg- 
ative regions on each side of the ring, looking like two negative 
rings inclined with respect to the positive ring. These two nega- 
tive regions cross the ring close to the minor and major axis and 
create truncatures, as in cADI. The intensity of these negative re- 
gions is significant: it reaches 50% of the ansae intensity for the 
reduced parallactic angle amplitude and 30% for the real (77.5°) 
amplitude. The angle between the negative and positive rings is 
directly driven by the separation criterion Ng and the number of 
images N used to build the reference PSF The smaller these pa- 
rameters, the closer the negative regions will appear with respect 
to the positive ring and the closer to the semi-major axis are the 
truncatures. 
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Fig. 6. Top: initial convolved disk and cADI-reduced disks for 
the real parallactic angle evolution (A9 = 77.5°) and for the re- 
duced amplitude (A8 = 23°). The color scales are identical and 
linear. Bottom: radial deprojected brightness profiles for the ini- 
tial convolved disk (dashed lines) and the cADI images (plain 
lines) in two directions (along the semi-major axis and 1 1 away 
from it). The color of the lines is the same as the color of the 
arrows in the top images. 




convolved disk at 1.1 arcsec 

cADI at 1 .1 arcsec 
cADI at 1.1 arcsec (reduced parang) 
convolved disk @ 1 .2 arcsec 

cADI at 1 .2 arcsec 
cADI at 1.2arcsec (reduced parang) 




40 60 
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Unlike cADI, rADI images are variable whether we observe 
the star before or after culmination when it passes the meridian. 
This creates a left/right and top/bottom asymmetry as shown in 
Figure [9] This is because the reference frame is built using a re- 
duced set of N images in rADI, here 20. The negative regions 
appear to be stronger on the right side when the star is observed 
before culmination, e.g. when the parallactic angle evolution rate 
is increasing. This effect is stronger for a small separation crite- 
rion. 

Discontinuity in the parallactic evolution pattern is another 
effect that create asymmetries detectable in rADI, as shown in 
figure [TO] It easily induces brightness asymmetries in the neg- 
ative regions up to 30%. The two negative regions are indeed 
symmetrical about the ring for the real observed parallactic an- 



Fig. 7. Top: deprojected disk. The two circles have a radius of 
1.1" and 1.2". Bottom: azimuthal brightness profiles for the ini- 
tial convolved disk, and the 2 cADI images along those two cir- 
cles of radius 1.1" and 1.2". An azimuth angle of 0° corresponds 
to the right in the top images. 



gle evolution, but this is no longer the case when a discontinuity 
is introduced at the beginning (discontinuity 1) or at the end (dis- 
continuity 3) of the sequence. When introduced in the middle of 
the sequence (discontinuity 2), the residual disk remains sym- 
metrical. 
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Fig. 9. Disks after rADI processing for the so-called "reduced" 
parallactic angle evolutions (AO = 23°). They reveal asymme- 
tries depending on observing conditions before/after meridian. 
The separation criterion Ng is 3FWHM for the two left images 
and 1.5 for the two images on the right. The color scales are 
identical and linear. 



Reduced par. angle 



Fig. 8. Residual disks after rADI processing for the observed 
parallactic angle evolution (AO = 77.5°) and the one with re- 
duced amplitude (AO = 23°). The separation criterion Ng is 3 
FWHM for the 2 images at the top and 1.5 for the 2 images at 
the bottom. The color scales are identical and linear. 



4.1 .3. Effects specific to LOCI 

In cADI or rADI, the way reference frames are built depends 
only on geometrical parameters. However, in LOCI, the linear 
combination also depends on the underlying speckle noise dis- 
tribution. As a result, a relevant analysis of disk images reduced 
in LOCI cannot be performed without noise in the data. 

For that purpos e, we used part of ou r L' data of HR4796 
from April, 6 2010 ( Lagrange et al.l2012bl) . in which we inserted 
a fake disk with a S/N of about 20cr to test a few realizations of 
the LOCI optimization algorithm. The data cube consists of 67 
images with 52° field rotation. Our goal here is not to make a 



comprehensive study testing all possible LOCI parameter sets 
and noise distributions, but to provide some hints on artificial 
effects occurring with LOCI. 

For the baseline LOCI reduction, we used the following pa- 
rameters: N s = 1.5 x FWHM, N A = 300, g = 1 and rfr = 
1 .4 x FWHM. To distinguish between effects due to the inher- 
ent noise of the data and effects due to the LOCI algorithm, we 
stored for each reduction the LOCI coefficients c k used to build 
the reference frames and applied them to a data cube containing 
the fake disk without noise. We show in FigureQT]the rea l LOCI 
result on the left and the reconstructed disk image if there was no 
noise in the data on the right, for each set of LOCI parameters. 

As in rADI, all sets of parameters show bright blobs in the 
ansae when the separation criterion Afo is small. The ring appears 
thinner elsewhere. The S/N is poorer close to the minor axis, 
therefore the noise can create distortion or intensity fluctuations 
in the ring that are not visible in the reconstructed image. 

No significant brightness asymmetry between the left and 
right ansae is detectable in the reduced LOCI images. Negative 
areas appear around and inside the ring, up to a level of 20% of 
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Fig. 10. Disks after rADI processing for the observed parallactic 
angle evolutions (A9 = 77.5°) and for the three different discon- 
tinuities introduced. The separation criterion N$ is 3FWHM for 
the two left images and 1.5 for the two images on the right. The 
color scales are identical and linear. 
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Fig. 11. Left: residual disks after LOCI processing of the real 
data cube of HR4796 taken on April, 7 2010 with a fake disk 
inserted at 90° from the real one. Different LOCI parameters are 
used and the coefficients are stored. Right: residual disk after 
applying those exact same LOCI coefficients to a data cube con- 
taining only the fake disk (without noise). The color scales are 
identical and linear for all images. The cADI and rADI reduc- 
tions are shown at the bottom for comparison. 



the flux in the ansae. The negative regions are more enhanced on 
the right side. 

The radial-to-azimuthal aspect ratio g appears to have an ef- 
fect on the intensity and width of the ring. Optimization areas 
more azimuthally extended (i.e. with lower g) result in less self- 
subtraction. TJmjs_akoJhe_casefor the set of conservative LOCI 
parameters dThalmann et al.l201 ll) . when Na = 10000, since this 
configuration means using the entire radial extension available in 
the image for each optimization area. The disk is 14% brighter 
in that case, but negative areas are also emphasized. 

The highest flux in LOCI restitution occurs when a low as- 
pect ratio g (here, 0. 1 which is the same as using a full annulus) 
is combined with many PSF cores (.^=1000). In this case, the 
disk intensity is almost similar to that obtained in cADI, except 
for the negative areas in the center of the ring. 

It also turns out from these LOCI-reduced images that it is 
difficult to predict with the LOCI algorithm. 



4.1 .4. ADI applicability with respect to the disk inclination 

So far, we have studied both the flux loss and the ADI artifacts 
with disks that have a fixed inclination. To answer the question to 
what extent in terms of disk inclination ADI can be applied, we 
performed additional cADI reductions of images of an inclined 
ring for several inclinations, as shown in Figure [12] Obviously, 
no clear answer can be provided since many additional aspects 
must be considered, such as the disk parameters to be retrieved 
and the amount of self-subtraction acceptable given the data S/N. 
However, we show here that for low-inclination, e.g. below 60°, 
the major and minor axis are critical regions affected by ADI 
because many low spatial frequencies are lost in those regions. 
Even for higher inclinations, the minor axis is always too much 
affected by self-subtraction to allow any information retrieval. 
Concerning the major axis, blobs at the tip of the ring ansae (de- 
scribed in more detail in Section I4.1.U appear below 70° there- 
fore great care must be taken not to misinterpret those features 
with physical processes such as dust grains blown by radiation 
pressure. These results suggest that disks inclined below 50° are 
inappropriate for ADI. 
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Fig. 12. HR4796-like annuli for different inclination (left) and 
the residual disk image after cADI for a total field rotation of 
50°. The color scale is identical and linear for all images. 
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Fig. 13. ADI reduction of a spiral arm seen pole-on, with and 
without noise from L data. The color scales are identical and 
linear for all images. 



4.2. ADI on spiral arms 

Recent disks disco veries included spiral a rms in circumstellar 
disks. For instance, Fukagaw a et alJ (120 06) performed reference 
star PSF s ubtracti o n to r eveal spiral arms around HD 142527, 
and iRameau et alJ d2012l) confirmed these features using both 
ADI and reference star subtraction. Other high-contrast imag- 
ing technique such as polarimetry were used for SAO206462 by 
iMutoet alJd2012h . 

Because spiral arms have a very specific geometry with az- 
imuthal spatial frequencies varying with the separation from the 
star, we carried out a specific investigation of the ADI impact on 
these features. 

An Archimedian spiral, seen pole-on, parametrized in polar 
coordinates by r = 10 x 9 with an initial width of 10 pixels was 
processed through cADI, rADI and LOCI with a total field ro- 
tation A9 = 37°. Results are shown in Figure [T3l To allow a 
relevant comparison of cADI and rADI images with LOCI im- 
ages, we added noise to the data. For that purpose we used our 
L data ( FWHM=3 .6px) of the M star HR 4796B taken on April, 
6 2010 dLagrange et al.ll2012bl) in which we inserted the spiral 
arm with an S/N of 14. 

The first conclusion is that arms are significantly thinner. The 
width reduction reaches 30% in cADI. 

As previously, we observe negative shadows on each side of 
the arm. These shadows can be deep and not necessarily sym- 
metrical about the arm. This is mostly visible in rADI without 
noise. When the curvature is smaller than the curvature of a cir- 
cle at the same location (e.g. for the innermost regions), the neg- 
ative parts on the outer edge of the spiral are enhanced, and vice 
versa. 



Great care should also be taken when interpreting the inten- 
sity distribution of a spiral arm. Although the initial spiral has 
a uniform flux, the LOCI and rADI reductions show a higher 
flux at short separation. The flux loss due to self-subtraction is 
indeed a trade-off between how large the amplitude of field ro- 
tation is and how small the azimuthal component of the spatial 
frequencies is. 

The case of a spiral arm inclined towards the line of sight 
was also investigated. It appears to be very similar to the case of 
an inclined ring. Some regions behave like the semi-major axis 
of an inclined ring, creating bright blobs when the amplitude 
of field rotation is low and some regions behave like the semi- 
minor axis of an inclined disk where most the disk flux is self- 
subtracted. 



5. Quantitative estimation of the biases 

We now provide a quantitative estimation of the change in some 
important disk observable properties from an astrophysical per- 
spective: the disk total flux, the outer slope of the SBD and the 
ring full width at half maximum (FWHM) measured radially 
along the semi-major axis. For that purpose, we focus on ring- 
shaped disks. We analyze the effect of the following variable 
parameters: the inclination of the ring, the radial exponents a out 
and a,„ that parametrize the volumetric dust density for three 
different ADI algorithms as detailed in Table Q] 

For each set of disk parameters and ADI algorithms, the re- 
ductions were performed with and without noise. For the case 
with noise, we again used part of our L' data of HR 4796A from 
April, 6 2010 (Lagrang e"etalll2012bl) . into which we inserted a 
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Table 1. Set of disks and ADI parameters used for the quantita- 
tive analysis. 



AO 


52° 


1 


60°;65 o ;70 o ;75 o ;80 o ;85° 


din 


2;3;5;7;10 


®out 


-2;-3;-5;-7;-10 




cADI 


ADI algorithm 


rADI N s = 1.5, N = 20, dr =1.4 
LOCI N s = 1.5, g = 1, N A = 300, dr = 1.4 



fake disk with a S/N of about 10<r. The data cube consists of 67 
images with 52° field rotation. 

The slope and FWHM measurements along the semi-major 
axis appear to be very sensitive to the noise, particularly for very 
narrow rings because we lack sufficient dynamic in the image. 
Therefore, to isolate the effect of ADI, we performed those two 
measurements on the images reduced without noise. For LOCI 
we computed the LOCI coefficients for data with noise and ap- 
plied them to a noiseless model of the disk. 

5.1. ADI impact on the total flux of the disk 

Following the general model detailed in Section I3.31 ADI be- 
haves as a high-pass filter, therefore the smoother the disk bright- 
ness distribution, the greater the disk flux loss. We illustrate this 
in Figure[14]in a number of quantitative examples coming from 
the reduction tests cases detailed in TableQ] 

1 . The amount of flux loss after ADI processing is very highly 
correlated to the inclination of the disk. The more pole-on the 
disk, the higher the flux loss, whatever the algorithm used. 
Figure [14]provides typical examples that can be used as ref- 
erences, with flux losses ranging from 100% to a few % for 
disks inclined from 60° to 85°. 

2. The inner and outer slopes of the dust density distribution 
also have a significant impact on the surface brightness of 
the processed disk, since the softer the slopes, the higher the 
flux loss. For the expected slopes of circumstellar disks, e.g. 
about -2 to -5 for a out and 2 to 5 for a,„, the flux loss varies 
between 50% and 80% , as visible in the middle and bottom 
chart of Figure [14] 

3. In cADI, flux losses are about 10% higher when the ref- 
erence frame is built using the mean rather than the me- 
dian. This conclusion contradicts the graphs of Figure [2] 
which were produced from noiseless data. The mean is in- 
deed much more sensitive to noise, which leads to more disk 
flux captured in the reference frame in the presence of noise. 

5.2. ADI impact on the outer slope of the surface brightness 
distribution 

Because the inner or outer slope of a debris disk is a significant 
parameter for understanding the physics of the disk, we studied 
the effect of ADI on the slope along the major axis of a ring. 
As already seen, the inner part of the ring is strongly affected by 
self-subtraction, therefore attempting to retrieve the dust density 
inner slope from ADI images is dangerous. For this reason, we 
focus here on the parameters that influence the outer slope of the 
SBD. 

The SBD outer slope is retrieved from a linear regression 
on a log-log plot of the SBD profile. Two elements have an im- 



pact on the slope of the final image: the convolution of the initial 
disk with the PSF that tends to make the slope of the profile 
softer, and ADI processing whose effect is non-linear. The rel- 
ative change in the SBD outer slope due to ADI reduction is 
presented in Figure [15] A negative change means a steeper pro- 
file. As the reduced LOCI image is no longer symmetrical about 
the semi-minor axis, we provide here the slope measurements on 
each side of the semi-major axis. We show here the impact of the 
parameters i and a ou , only, as the inner slope of the dust density 
a ln has no influence on the SBD outer slope. 

1. in cADI (median-combined reference frame), the inclination 
and initial value of a out has little impact on the slope. The 
slope is slightly softer for values of a out close to -2. 

2. in LOCI, the effect is variable whether we consider the East 
or West side of the ring. The trend seems however to follow 
that of cADI. 

3. in rADI, the SBD slope is clearly steeper, by -20 to -5% for 
the explored range of inclinations and a out . 

On the whole the effect is all the more pronounced as the disk is 
spatially extended with low inclination and soft SBD slope. 

5.3. ADI impact on the FWHM of the ring 

In almost all simulation cases, the measured FWHM is slightly 
decreased after ADI reduction. 

1. First, the inclination of the disk can significantly bias the 
measured FWHM. The more face-on the disk is seen, the 
higher the impact on the FWHM compared to the initial con- 
volved disk (see Figure[l6]top chart). 

2. The measured FWHM of a ring is also highly dependent on 
the dust distribution inner slope of the initial disk. Indeed, for 
very soft radial density exponents a out , the initial FWHM of 
the ring is large. ADI then creates negative areas inside the 
ring, consequently the inner half-width is much decreased. 
This effect is visible in the middle chart of Figure [TBI 

3. Last, the dust density outer slope has a fairly low impact on 
the disk FWHM. In rADI, the change in FWHM is constant 
over the explored range of a our and the change decreases 
from -38% to -22% (without noise) when a ou , approaches 
(Figure[T_6] bottom chart). 

5.4. Sensitivity analysis relative to the LOCI parameters 

The study presented above compares the reduction of various 
initial disks but with a single set of rADI and LOCI parameters. 
To complete this study, we now show for a given disk morphol- 
ogy, namely a ln = 5, a out - -2.5 and i = 75°, a comparison of 
different LOCI parameters. 

We have already shown that the separation criteria Ng is the 
main parameter that influences the amount of self-subtraction 
of a disk. The size of the LOCI optimization area, expressed in 
number of PSF cores Na, is the second contributor. The larger 
this area, the lower the flux loss, as shown in the top chart of 
Figure [T7] For a large domain, the chances that the RMS opti- 
mization matches the disk pattern and not the speckle pattern is 
indeed lower. 

Now, for a constant size of the optimization area (here Na 
is kept to 300), the aspect ratio is also influencing the amount 
of self-subtraction, as shown in the bottom chart of Figure [17] 
For the morphology studied, the larger the radial extension of 
this area, the less self-subtraction occurs. We interpret this fea- 
ture by the fact that the relative area contaminated by the disk is 



10 



J. Milli et al.: Impact of angular differential imaging on circumstellar disk images 




60 65 70 75 80 

Inclination i in ' 




2 4 6 8 10 

Inner slope of the dust density a in 




Fig. 14. Flux loss on a disk as a function of the disk inclination i (top chart), the dust density inner slope ff,„(middle chart) and the 
dust density outer slope a oM (bottom chart). The flux loss is the ratio of the total flux of the reduced disk by the total flux of the 
initial disk after convolution by a synthetic PSF. For the top chart, a,„ and a out are kept constant to 3 and -5, for the middle chart i 
and a out are kept constant to 75° and -5, and for the bottom chart i and a,„ are kept constant to 75° and 3. 



lower for the optimization zones characterized by low values of 
g. This conclusion is obviously valid for the geometry of disks 
considered here and should be revisited for a disk different from 
an inclined ring. 

5.5. Position angle and offset of the disk 

The position angle (PA) and offset of a disk are very important 
parameters for describing the morphology of a disk. Within ADI, 
two different reasons can alter the PA and offset as measured on 
the reduced disk image: 

- star centering 

- artifacts resulting from disk self-subtraction. 

Star centering is an important problem in ADI because image 
rotations or derotations are frequently performed, and one does 
not always accurately know the star center for coronographic or 
saturated images. The PA of a disk is less sensitive a parame- 
ter than the PA of a planet, however, as Lagrange et afj (1201 2al) 
showed for the f3 Pictor is disk and its planeta ry companion. In 
the case of HR4796A (Lagrange et al.ll2012bl) . the disk PA un- 
certainty due to the imperfect knowledge of the NaCo PSF center 
reaches 0.2°, which is the same order of magnitude as the uncer- 
tainty associated to the PA measu rement technique (th e maxi- 
itium regional merit m ethod, as in iBuenzli et al.l (l2010h and in 
Thalmann et al.l (1201 ll) ). Concerning the the uncertainty on the 
disk cente r due to the uncertainty on the PSF center, we also 
showed in lLagrange et alJ (1201 2al) that it is the same order of 
magnitude as the uncertainty on the star center, 0.2 px in that 
case. 



Then we investigated the influence of artifacts resulting from 
the disk self- subtractio n on the PA and offset measurement. In 
Lagrange et al.l (1201 2bl) . we reduced simulated disk images of 
HR4796A in cADI, rADI and LOCI and compared the PA and 
offset of the reduced disk image to their inital value, as measured 
with the regional merit method. Two different initial disk mod- 
els were used, called HR4796SD and HR4796blowoutSD, and 
77.5° field rotation was used in the simulation. In both cases, the 
error on the disk PA due to ADI was below 0.1° or 0.1%. The 
measured bias on the offset of the disk was within the measure- 
ment uncertainty of the regional merit method. 

We emphasize that some specific configurations might nev- 
ertheless bias the PA measurement. We see indeed in Figure [8] 
that rADI reduction with 23° field rotation tends to artificially 
create a discontinuity between the right and left sides of the ring 
that can bias the PA measurement if made separately for each 
side. 

5.6. Overestimation of planet detection limits 

It is well known that dust sur rounding stars is a p ossible source 
of confusion with exoplanets (Roberge et al.ll2012h and can lead 
to false-positive detections. Even without detection, the presence 
of a disk can also lead to an overestimation of planet detection 
limits. 

The detection limits are usually computed by first calculating 
the flux loss for a fake planet injected in the data and then apply- 
ing this correction factor to an estimation of the 5cr noise level in 
the reduced ADI image. Great care should be taken in this case 
because the region inside the annulus and around it has a neg- 
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Fig. 15. Bias on the measured outer slope of the disk SBD as a function of the disk inclination i (top chart) and the dust density 
outer slope a om (bottom chart). For the top chart, <?,„ and a out are kept constant to 3 and -5 and for the bottom chart, i and a ln are 
kept constant to 75° and 3. 
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ative flux, as demonstrated previously and illustrated in Figures 
[8] [9] and [10] Therefore, any planet located in this region and that 
has a flux comparable to the disk flux will undergo severe flux 
losses from the disk itself in addition to its own self-subtraction. 
The position and amplitude of these negative areas are difficult 
to predict, therefore planet detection limits can be easily overes- 
timated. 



To give a quantitative example of this effect, we used the 
model of HR4796A disk and injected two fake faint planets at 
two critical locations where much flux loss from the disk self- 
subtraction is expected (see Figure [T8l>. The planet fluxes were 
scaled relative to the disk brightness to correspond to a \AMj UD 
planet according to the evolutionary model of Baraffe et al. 
(120031 or a flux ratio between the planet and the disk in the 
ansae of around 1:2. Then ADI reduction was performed in 
cADI and rADI for two amplitudes of field rotation AO = 77.5° 
and 23°. The rADI parameters are a separation criterion Ng = 1.5 
FWHM, a number of images N - 2Q and an annulus width 
dr — 1 .4 FWHM. No noise was added to the data and we con- 
sidered a baseline scenario of two planets alone without disk. 
We compare the results in terms of planet flux loss in Table [2] 
The presence of the disk significantly impacts the detectabil- 
ity of both planets because the reference frame captures much 
flux from the disk at the position of the two planets. Although 
the planet self-subtraction is always well below 10% even for a 
small field rotation, the presence of the disk very often leads to a 
100% flux loss, meaning that the planet is no longer detectable. 



Fig. 18. Position and brightness of the two fake planets relative 
to the disk model of HR 4796A. The color scale is linear. 

6. Improvements to ADI methods: masking, 
iterating, and damping techniques 

Several ADI techniques more adapted for extended geometries 
are reviewed in this section. The biases they involve on disk ob- 
servables are reported in Table [3] and are commented on in the 
following subsections. 

6.1. Iterations 

The main challenge of applying ADI on disks is to separate the 
disk from the star to build reference frames. The very first so- 
lution proposed to tackle this problem is an iterating process, 
applied to reveal the /3 Pictoris disk from HST/STIS images 
(Heap et al. 2000). It was reapplied more recently for the disk of 
HD 207129 on HST/ACS images dKrist et al.ll2010l) . For those 
two space-based data sets, only two images at two different roll 
angles were used. However, one can apply this principle for a 
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Fig. 16. Bias on the FWHM of the ring (measured along the semi-major axis) as a function of the disk inclination i (top chart), the 
dust density inner slope ff,„(middle chart) and the dust density outer slope a out (bottom chart). For the top chart, a m and a out are 
kept constant to 3 and -5, for the middle chart i and a om are kept constant to 75° and -5, and for the bottom chart i and a m are kept 
constant to 75° and 3. 



Table 2. Flux loss from two fake planets after ADI reduction. The flux loss is computed by aperture photometry with a radius of 
0.75 x FWHM. mcADI is a technique introduced and commented on Section lo*T2l 





cADI 


mcADI 


Ad = 11.5° 
rADI Ns = 1.5 


rADI N s = 3 


cADI 


mcADI 


A<9 = 23° 
rADI Ng = 1.5 


rADI N 6 = 3 


Planet 1 


without disk 
with disk 


0.2% 
27.1% 


0.3% 
4.3% 


1.4% 
100% 


0.2% 
23.5% 


8.5% 
100% 


NA 
NA 


2.7% 
100% 


2.6% 
99.7% 


Planet 2 


without disk 
with disk 


0.1% 
4.4% 


0.3% 
1.8% 


1.4% 
100% 


1.2% 
100% 


3.1% 
63.0% 


12.9% 
21.1% 


2.7% 
45.0% 


0.7% 
100% 



Table 3. Summary of the disk biases with standard and disk-specific ADI reduction strategies. Flux losses and FWHM changes were 
measured on data with and without speckle noise respectively, for an initial disk parametrized with i = 75°, a m = 5 and a ou , - -2.5. 
The rADI and LOCI parameters are Ng = 1.5, Na = 300, g = 1, except for conservative LOCI, which has Ng = 2, Na — 1000 and 
S = 0.5 





cADI 


mcADI 


cADI 
iter. 1 


cADI 
iter. 2 


rADI 


LOCI 


mLOCI 


cons. 
LOCI 


LOCI 
nnls 


dLOCI 
A = 1 


dLOCI 
A = 10 


Flux loss (%) 


60.7 


21.6 


53.4 


51.6 


69.1 


71.8 


56.3 


53.8 


66.7 


64.3 


49.0 


FWHM 

change (%) 


-21.7 


-3.3 


-0.3 


-2.4 


-32.7 


-24.9 


-21.9 


-13.6 


-31.4 


-29.5 


-21.2 



sequence of ground-based images take n in pupil-tracking mode . 
It is described for a single iteration in Lagrang e et al.l d2012al) . 
A first estimate of the reference PSF image is made by median- 
combining all images and subtracting them from the images as 
in cADI. The residual image after derotation and collapse is the 
first estimate of the disk. In a second step, the contamination of 



the reference frame by the disk is computed using the estimated 
disk image of step 1 . This contamination is then subtracted from 
the reference frame of step 1. PSF subtraction, derotation and 
collapse then yield the residual disk image of step 2. This pro- 
cess is then repeated iteratively. This equals applying sequen- 
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Fig. 17. Flux loss on a disk as a function of the size of the LOCI optimization region Na (top chart) and its radial-to-azimuthal 
aspect ratio g (bottom chart). An initial disk parametrized with i = 75°, ai„ = 5 and a out = -2.5 is used here, with g kept constant 
to 1 . for the top chart and Na kept constant to 300. for the bottom chart. Flux losses in c ADI and rADI (Ng = 1 .5) are shown here 
for reference. 



tially cADI after correcting the reference frame from the residual 
image of the previous iteration. 

Simulations on fake disks show that this process helps in- 
deed to recover part of the disk flux lost by the effect of self- 
subtraction. More importantly, it can reduce the amplitude of 
the negative areas in the reduced image if one uses the a priori 
information that the estimate of the disk image must be posi- 
tive. However, solutions converge toward a final disk image that 
is still less extended than the initial disk since information is 
lost in the areas where negative flux appears after the first guess. 
This process was applie d to two real astrophysi cs ground-based 
observations: 8 Pictoris (Lagrang e et al.ll2012al) and HR4796A 
(Lagrange et alJl2012bl) . Because the amplitude of field rotation 
was large, the gain was not significant and convergence was al- 
ready reached after two iterations. 

It was described here for cADI, but it can also be used with 
rADI, the difference being that there is one reference frame for 
each image of the data cube. 

The typical flux recovery that can be expected for a disk in- 
clined at 75° with a m = -2.5 and a out = -5 with one and two 
iterations are given in Table [3] It shows that around 10% more 
flux is preserved after the first two iterations. Close to the star, 
we emphasize, however, that it also amplifies residual structures 
that are not part of the disk, such as bright speckles or residu- 
als of Airy rings. Far from the star, this techniques turns out to 
be efficient to preserve the initial disk surface brightness, since 
the measured FWHM along the semi-major axis is only slightly 
affected with changes smaller than 3%. 



6.2. Masking the disk 

One of the drawbacks of iterating techniques is that if a bright 
speckle is present in the reduced disk image after the first it- 
eration, it will be attributed to the disk extended emission. 
Therefore, it will stay there in the following iterations. One way 
to circumvent this difficulty is the following: instead of using the 
first iteration result as an a priori of the disk emission to subtract 
its contribution in the reference frame, one can also build an a 
priori binary mask of the disk and take into account only pixels 
outside this mask when building the reference frame. This obvi- 
ously reduces the S/N, especially for very close regions, but the 
gain in outer regions, e.g. in the ring ansae, can be significant and 
negative areas are also li mited. An example of ma sking in cADI 
is shown for HR 4796 A in Lagra nge etall (1201 2bl) . It shows that 
the ansae are 10% brighter with masking. 

Many aspects make this algorithm particularly attractive for 
disk reduction: 

1. The flux of the disk is well preserved for a well-defined 
mask (an additional 40% of disk flux is recovered with re- 
spect to standard cADI, as shown in Table [3). Moreover, be- 
cause negative regions are limited, the risk to overestimate 
planet detection limits is lower. This is shown in Table[2j for 
a sufficient amplitude of field rotation, the flux loss is limited 
to 4.1% for the closest planet (planet 1), whereas it reaches 
27.1%incAD0 

3 We can also notice that flux losses are slightly higher in mcADI 
than in cADI without disk (0.3% vs 0.2% for planet 1 and 0.3% vs 
0.1% for planet 2). This is because when masking is applied, there are 
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2. The way the images are combined to build the reference 
frame is deterministic and independent of the noise distri- 
bution 

3. All available images are used to build the reference frame, 
therefore fewer artifacts as in rADI are expected. 

4. Using a disk binary mask is more relevant than using a sepa- 
ration criterion Ns since this means adapting the criterion to 
the spatial extent of the disk. It requires that the disk is not 
too extended, however. 

Masking the disk can also be applied with the LOCI re- 
duction. Taking into account a bin ary weighting fun ction in the 
LOCI algorithm is described in IPuevo et alJ (1201 21) for point- 
source detection. It can be used for extended objects as long as 
the size of the optimization area after withdrawing masked pix- 
els Na is still large enough for the inverse problem to be well- 
conditioned (Na > Nr, where Nr is the number of available 
reference images). As displayed in Table [3] the biases on disk 
observables have similar orders of magnitude as for a standard 
LOCI algorithm when conservative parameters are used (Ns - 2, 
Na = 1000 and g = 0.5). 

6.3. Damped LOCI 

Damped LOCI is an algorithm that forces the coefficients used to 
construct the reference frame to be positive (referred to as LOCI 
nnls in Table [3] for non-negative least squares) and uses a mod- 
ified cost function with respect to standard LOCI to compute 
them (IPuevo et al.ll2012l) . It seeks to minimize the residual flux 
in the optimization region while maximizing the flux in the sub- 
traction region. The maximization constraint is set by introduc- 
ing an additional penalty term in the cost function. The relative 
importance of this term depends on the Lagrangian multiplier 
factor A. Therefore the final reduced image will be a trade-off 
between minimization of the speckle noise and preservation of 
the disk flux. It has only been applied to point-like sources so far 
and we propose to investigate its effect on extended objects. 

Figure [19] shows that the constraint on the positivity of the 
LOCI coefficients has a strong effect on the preservation of the 
disk, with an 8% gain on the recovered disk flux. It also performs 
better than mLOCI. The impact of increasing the damping coef- 
ficient from to 1 on the recovered disk flux is linear but rela- 
tively small (3% additional gain). An additional increase in A up 
to 10 enables one to gain even more flux at the expense of the 
speckle noise attenuation. On the measured FWHM, changes in- 
duced by dLOCI with A < 3 are slightly higher than for standard 
LOCI and reach the same level as cADI for A = 10. A qualita- 
tive comparison of disks reduced with standard LOCI, mLOCI 
and dLOCI for A = 1 and A = 10 is provided in Figure [20] As 
explained, this shows that more disk flux is preserved in mLOCI 
and dLOCI. However, we also clearly see that speckles pinned 
on the diffraction rings are less well subtracted in those cases. 
Bright blobs described in Section |4. 1 . II appear in the two ansae 
for high damping coefficient (dLOCI, A = 10) or in mLOCI. 



7. Conclusions 

We discussed the effects of ADI applied to extended circum- 
stellar disks. A theoretical study demonstrated that the field of 



fewer frames available for the reference frame construction. Therefore, 
the relative number of frames with a small amplitude of field rotation is 
higher in mcADI than in cADI. 




Fig. 20. Right: comparison between standard, masked and 
damped LOCI (A = 1 and A = 10) for the same N s = 1.5, g = 1 
and Na = 300 on a disk parametrized by a/,, = 5, a out = -2.5 
and i = 75°. Right: corresponding disks after applying the same 
LOCI coefficients to a data cube containing only the fake disk 
(without noise). 

application of this imaging technique depends on two parame- 
ters: the available amplitude of field rotation and the geomet- 
rical extension of the object in azimuth. This analysis gives a 
clear indication of the required amplitude of field rotation to re- 
duce a given disk depending on its morphology, and conversely 
it gives an indication of the disk minimum inclination for a given 
field rotation. With reasonable observing time, objects inclined 
less than 50° are out of the scope of ADI. For all others, self- 
subtraction is a problem that observers must bear in mind when 
reducing the data. It affects the whole geometry of the disk by 
reducing the total disk flux, but also introduces local effects dif- 
ficult to distinguish from speckle noise patterns. A typical effect 
already observed for some circumstellar disks reduced with ADI 
is the enhancement of the ansae for an inclined ring morphology, 
creating bright blobs along the disk major axis that can be easily 
misinterpreted. Brightness asymmetries from non-uniform par- 
allactic angle evolution are other striking effects that arise from 
ADI, especially rADI in which a limited number of images is 
used to build the reference frames. 

Measurable astrophysical quantities used to characterize 
disks such as inner/outer slopes of the SBD or ring FWHM are 
biased and their estimation must be calibrated. We here provided 
an insight into how these parameters are affected, depending on 
the initial disk morphology and the reduction parameters. 

To evaluate the efficiency of the explored ADI algorithms ap- 
plied on disks, two key points must be considered: the amount 
of disk self-subtraction and how well the residual speckles are 
attenuated. For the first point, masking as used in mcADI is our 
preferred data reduction strategy, since it replaces the classical 
separation criterion Ns appropriate for point sources with a bi- 
nary mask adapted to the expected extension of the object, to 
efficiently limit self-subtraction. The second point is better ad- 
dressed with optimization algorithms such as LOCI, but flux 
losses are much more severe. Our study shows that using larger 
optimization areas, with an azimuthal extent greater than the ra- 
dial extent, helps to preserve the disk flux. More importantly 
constraining the LOCI coefficients to be positive proves to be 
very effective in this goal. 
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Fig. 19. Impact of the damping coefficient A on the disk observables. The value of the parameters N$, g and Na are kept constant 
to 1 .0, 1 .0 and 300, respectively. The flux change in cADI, rADI and standard LOCI is displayed for reference, the x-axis is on a 
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Appendix A: Deriving the flux loss criteria in cADI 

We derive the expression of the residual flux given in EquationQ] 
which is valid for a disk seen edge-on, with an intensity vertical 
profile parametrized by the expression I(z) = e -(f)' . Such a disk 
is plotted in Figure lATTl The vertical scale height h is in this case, 
constant with the radius, but the following derivation still holds 
if h varies with the separation to the star, for example if there is 
a flaring characterized by a flaring index /3 such that h = h(x) = 

Similarly, the intensity in the disk mid-plane (z - 0) is 
constant here and equal to 1 , but the derivation still holds if the 
intensity varies with the separation x. 

We consider a point A located at a separation x and a height 
Z. It makes an angle i = arctan with the disk mid-plane. The 
vertical profile at this separation as a function of i is displayed 
in Figure IA.3I The two blue arrows represent the amplitude of 
field rotation A9/2 on each side of A. In cADI, the flux of the 
reference frame at the position of A is equal to the median of 
the curve between i - A9/2 and i + A9/2, sampled according to 
the parallactic angle distribution. We assumed that this sampling 
is uniformly distributed. To compute the median, we divide this 
range into four regions, according to the values of the intensity: 

- region 1, of length 2i, the farthest on the right characterized 
by the lowest intensity values 

- regions 2 and 3, symmetrical about the mid-plane, of length 
A0/2-2z 

- region 4, of length 2i centered on the mid-plane and consist- 
ing of the highest intensity values. 

This decomposition into four regions of increasing intensity 
allows us to directly conclude that the median over the four re- 
gions is also equal to the median over region 2 only, which is 
equal to the intensity at point B for Ad/ 4. 

The disk flux / after subtraction of the reference frame can 
then be expressed as 

7(r, Q - 7(r, A9/4) 
; ~ Kr,i) 

_ »r(-(£) y )-»p(-(ffi) y ) 

= l - eXP [ Aa? j 



Initial disk 




cADI (median) 




Fig. A.l. Theoretical edge-on disk. 



An experimental verification of this law was performed and 
is shown in Figure [A3] The disk flux / in the mid-plane (i = 0) 
after cADI reduction, expressed in percentage for various disks 
and field rotations, is plotted against the criterion k c . All curves 
match well, which confirms the derived theoretical expression. A 
similar graph can be made for different vertical profiles instead 
of mid-plane profiles. 
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Fig. A.3. Experimental verification of EquationQ~|on simulated disks (with different values of the vertical exponent y and the vertical 
height h) for various field rotations AO. For k c = 3, 95% of the disk flux f is preserved. 



